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NOTATION 


"b wing span 

b^(ri) quadrature coefficient 

c local chord 

rate of change of hinge moment coefficient with flap deflection angle 
measured at a = 0^ 

rate of change of rolling moment coefficient with flap deflection 
angle measured at a = 0° 

Cl^ rate of change of lift coefficient with flap deflection angle mea- 

sured at a = 0° 

Ac^ change in section lift coefficient due to flap deflection 

AC^ change in overall lift coefficient due to flap deflection 

AC^ differential pressure coefficient 

E(m) complete elliptic integral of the second kind with parameter m 

h function that need not be precisely defined for the purpose of the 

section in which it is used 

h^(C) chordwise pressure mode 

J number of stations used in performing a numerical integration 

K(m) complete elliptic integral of the first kind with parameter m 

M Mach number; also an integer representing the number of integration 

points used in some particular numerical quadrature 

N number of chordwise control points 

x^(fi) position of the hinge line or its extension measured relative to the 
50-percent chord line and made dimensionless by the local semi- 
chord; -1 < Xq < +1 

a streamwise slope of the wing mean camber surface 

e /I - 

(S difference or an error or a control surface deflection angle 

6q constant used to control the overall accuracy of the chordwise 

integrations 


iii 



e quantity meant to te a constant or to shrink to zero in some process 

ni,ri2 spanw’ise positions of the flap side edges |ri2l ^ 1^1 I Q-^d r \2 > 0 

s\/eep angle of the control surface hinge line 

coordinates made dimensionless hy the semispan, see figure 1 (primes 
are used when the coordinate is an integration variable) 

symbol implying proportionality 

Subscripts 

c,f control surface or flap 

j particular spanwise station 

£e leading edge 

te trailing edge 


iv 


STEADY, SUBSONIC, LIFTING SURFACE THEORY FOR WINGS WITH 
SWEPT, PARTIAL SPAN, TRAILING EDGE CONTROL SURFACES 
Richard T. Medan 
Ames Research Center 

SUMMARY 


A method for computing the lifting pressure distribution on a wing with 
partial span, swept control surfaces is presented. This method is valid within 
the framework of linearized, steady, potential flow theory and consists of 
using conventional lifting surface theory in conjunction with a flap pressure 
mode. The cause of a numerical instability that can occur during the quadra- 
ture of the flap pressure mode is discussed and an efficient technique to elim- 
inate the instability is derived. This technique is valid for both the flap 
pressure mode and regular pressure modes and could be used to improve existing 
lifting surface methods. Examples of the use of the flap pressure mode and 
comparisons among this method, other theoretical methods, and experiments are 
given. Discrepancies with experiment are indicated and candidate caioses are 
presented. It is concluded that the method can lead to an efficient and accu- 
rate solution of the mathematical problem when a partial span, trailing edge 
flap is involved. 


INTRODUCTION 


Linearized, subsonic lifting surface theory for predicting lifting pres- 
sure distributions and the associated loads and moments on wings having arbi- 
trary planform and smooth surface slope distributions has been well developed. 
This is not the case for wings having surface slope discontinuities, however. 
Various methods for solving the control surface problem have been devised, but 
all seem to possess inadequacies that result from not accurately portraying the 
pressure distribution near the corner and at the side edges of the control sur- 
face. Landahl (ref. l) has determined the form and strength of the singular 
parts of the pressure distribution. This report will show how his results, 
suitably modified for a swept flap and specialized to the steady case, were 
incorporated into an existing lifting siorface theory. These results can also 
be easily incorporated into lifting surface methods based on the theories of 
Multhopp (ref. 2), Truckenbrodt (ref. 3), and Hsu and Weatherill (ref. k). 
Furthermore, the methods used herein indicate how the accuracy and efficiency 
of existing, conventional lifting surface methods can be improved. 


BASIC APPROACH 


The fundamental equation of linearized, steady lifting surface theory is 

^te^n’) 

f zAC (V ,n') 

2 K[5-C',B(n-n’)]dC'dn' (l) 

(n-n* )^+z^ 

where 

K[g-g* , e(n-n* )3 = 1 + ' (2) 

/(5-C )2 + 62(n_ri. )2 + g2z2 
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Figure 1.- Definition of coordinate 
system. 


(see fig. l). This equation is the 
mathematical expression of the fact 
that the wing is represented by a dis- 
tribution of semiinfinite line doub- 
lets with strengths proportional to 
ACp and with their end points placed 
on the wing and at downstream infinity 
(see ref. 5 for a derivation). In the 
analysis problem, a(S 5 n) is given and 
ACp(?’,n’) must be determined. The 
limit and differentiation processes 
can be interchanged with the integra- 
tion if special procedures are 
observed when the integrations are 
performed. Thus the integral equation 
is usually written as 




i(5 




ACp(5’, n’) 


K[5-5' ,6(n-n’ )]d5’ dn' (3) 


(n-n’)2 

(The cross on the integral sign symbolically designates these procedures.) 




The method for solving this integral equation for a wing of arbitrary 
planform is to first express the \mknown ACp(C’,n') in terms of known func- 
tions with unknown coefficients. The functions generally used are obtained 
from a consideration of two-dimensional airfoil and lifting line theories. A 
convenient expansion, based on the theories of Multhopp and Truckenbrodt , is 
that used by Wagner (ref. 6): 


ACp(?', n') 


N 


n=l 


(U) 
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The prescribed functions h^ are taJ^en from two-dimensional airfoil theory and 
are comprised of the familiar cot sin (j), sin 2(j), . • •, terms. The fn 

terms are still unknown, but can be prescribed in terms of known functions 
(which originated from lifting line theory) with unknown coefficients by the 
use of Multhopp’s interpolation formula: 


where 
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(-l)“‘^^sin e sin[(M+l)e'] 
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(M+l)(n’-n^) 
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( 6 ) 


Note that is not singular (fig. 2). Substituting equations (^) to (6) 

into (3) leads to 

N M 

( 7 ) 


where 


in) = Z) S 


n) = ^ 


nra 


n=l m=l 



te ACp iV i n' ) 

-iiS K[c-5',6(n-n')]d5' dn' (8) 




n 


m 


(9) 


Figure 2.- Multhopp interpolation 
function of equation (6) for M = 5 
and m = 2. 


The final step in the basic solution 
is to determine the unknowns , fnm. 'by 
enforcing equation (j) at a set of 
selected control points on the wing 
(collocation). 

The ownin’ s are called downwash 
modes (or regular downwash modes). 
Because of the particular choices for 
the pressure modes given in equation 
(9), the resulting downwash modes are 
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smooth functions (except at spanwise stations where c(ti’) is not continuous or 
differentiable), A moderate number of such downwash modes is inadequate to 
represent the discontinuous a(^ 5 n) which would occur for a wing having a 
deflected control surface. 


The most obvious method oT hand- 
ling a discontinuous a(Cjn) is to use 
the results of two-dimensional airfoil 
theory for a section with a flap (ref, 
T) and lifting line theory for a dis- 
continuous downwash distribution (ref. 
8) to obtain an additional (flap) 
pressure mode in the form of equation 
(i^). It was concluded in reference 9, 
however, that such a pressure mode 
cannot be made to reliably represent 
the partial span flap, and that pres- 
sure modes developed from three- 
dimensional theory are needed. 

Landahl (ref. l) developed such a 
pressure mode (fig. 3) for an unswept 
flap. He also showed that the pres- 
sure mode produces a downwash mode 
with exactly the required discontinu- 
ity at the leading and side edges of 
the flap, except possibly near the trailing edge. It has been determined that 
the pressure mode will produce the correct discontinuity at the trailing edge 
also. However, only the strength of the singular part of the pressure was 
determined by Landahl and a residual part remains to be determined. It is 
apparent then that the pressure mode need not contain any unknown elements nor 
be involved in a collocation procedure. The proper manner in which to employ 
Landahl*s result is quite simple. Let the integral equation (eq. (3)) be 
symbolically denoted by 

a = AC * K 
P 

Now let AC^ = AC-n + AC . - AC then 
P P pf pf * 

a = AC ^ * K + (AC - AC ^) * K 
pf p pf 

= a_ + (AC - AC ^) * K 
f p pf 

or 

ot - = (AC - AC -) * K 

f p pf 



Figure 3.- Distribution of the flap 
mode of pressure on an aspect 
ratio 2 rectangular wing with a 
single flap; hinge line at 
5 = 0.65 and edges at n = 0 
and T) = 0.5. 


The term 


can be calculated since 


AC 




is a known function and since 


is given, a-a^* (which could be termed a residual downwash) is completely known. 

The residual pressure distribution, AC -AC is unknown, but since a-a^ is 

P pf f 

continuous (although it may not be differentiable at a flap edge), AC^-AC^^ 
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should he sufficiently smooth that it may he expanded in the regular pressure 
modes. Consequently, the residual pressure may he determined with ,an existing 
lifting surface program, and the final pressure distribution is given merely hy 
adding the flap pressure mode to the residual pressures. This is basically the 
method followed here. The most difficult part was the numerical integration 
proced-ure used to calculate the flap downwash distribution efficiently and 
accurately. 

The procedures used are discussed in the next section in sufficient gen- 
erality that they may be used to compute the downwash caused by any pressure 
mode. The actual form of the flap pressure mode used is discussed in appendix 

A. 


BASIC INTEGRATION PROCEDURE 


It was necessary to adopt an optimized quadrature procedure for the flap 
pressure mode for the following reasons: 

(1) Evaluating the flap pressure mode at a single, given point is more 
difficult than for the regular pressure modes. 

(2) The pressure is logarithmically singular at the hinge line. This 
necessitates an increased number of chordwise quadrature points. 

(3) More spanwise integration points are needed than for regular pressxire 
modes and, due to the second-order singularity in the spanwise integral, the 
number of chordwise integration points must be increased simultaneously. 
Therefore, if the mmiber of spanwise integration points were doubled, the total 
number of integration points should be more than doubled (approximately 
tripled) . 

(U) The flap pressure is not in a separated form since it was derived from 
a three-dimensional analysis. Therefore, the pressure at each point must be 
calculated directly and cannot be represented as a simple product of previously 
calculated and stored values, that is, the matrix of the flap pressure mode 
values is not the outer product of two vectors as it is for the regular 
pressure modes. Consequently, the time required to perform a calculation for 
the flap pressure mode is much greater than for a regular pressure mode and 
therefore more efficient quadrature procedures are desirable. 

To best explain these procedures , consider first the spanwise integration. 
The flap downwash mode or any regular downwash mode can be manipulated into the 
following form: 


T >1 n, n’) 

^ 

^ (n-n’)' 


( 10 ) 


where Gf is the result of a chordwise numerical integration that has already 
been performed and may contain some error. The integration indicated by equa- 
tion (10) may be performed by various methods; the one used in this report is 
from Multhopp (ref. 2 ; see appendix B also), but the general procedure for 


5 


obtaining an efficient method would be similar for other choices of the span- 
wise integration technique. From reference 2, the following approximation is 
obtained for equation (lO): 

J 

n) = I t'.(n)G (c, n, n.) (ll) 

0=1 ^ 

where 


and 


n. = cos 9. 
J J 


9. = 


JJL 


j J + 1 


If n is restricted to 


n = cos 


VTT 


J + 1 


where v is some integer, then the expression for b (n) is 

J 


r 


bj(.) = ( 


j+i _ , 

4 sin 0 ’ ^ ^ 

J 

-sin 0 

^ > 1 V- j 1 odd 


(j+l)(n-ilj) 


0 , otherwise 


Vw 


( 12 ) 


Nonzero values of 4)j(n) for various choices of J are given in figure 4. 
Note that the coefficients become very large fop the integration stations 
nearest n. 


Consider the application of equations (ll) and (l2) for p = 0; 


af(C.O) , I , 

J -1 >9 »2 


sin 9 

f G-(0,C,n.) 

(nj)^ 


Recall that Gf is the result of a chordwise integration subject to some 
error and assume, for illustration, that it is calculated exactly at all nj 
except j = v-1, where the magnitude of the error is 6Gf. 


Then the magnitude of the error in computing is 


6a 


f 


-sin 0 

6G 

(j+l)(n^__^) 


which, for moderate or large J, is approximately 


6 



«“f = - ^ 6Gj, (13) 


-bj(T7) j 



Figure 4.- Quadrature coefficients for 
the downwash integral n = 0. 


Thus, by increasing J, the value 
obtained for af(C,0) does not con- 
verge to the correct ansver, but, on 
the contrary, diverges from it. 

A good example of this type of 
instability is given by Rowe (ref. 

10 ). He calculated the lift-curve 
slope of a square wing using various 
numbers of chordwise and spanwise 
integration points. When he used U 
chordwise points and increased the 
number of spanwise points beyond 12, 
his results diverged, whereas with 13 
chordwise points the results were 
stable when the number of spanwise 
points was increased to l6. If Rowe 
had continued increasing the number 
of spanwise integration stations 
without increasing the number of 
chordwise points, however, his 
results would have again become 
unstable. 

If follows from the discussion 
relating to equation (13) that, to 
guarantee convergence, the chordwise 
integration that determined Gf 
should use enough points so that 


6G 


f 


e(c, 0, n^_^) 

(J + 1)^ 


(14) 


where P is a constant greater than 1 (a value of 2 is used here). 


In effect, equation (l4) is satisfied typically by selecting the largest 
value of J likely to be used (or some equivalent to J if a quadrature rule 
other than that given by eqs. (ll) to (l2) is employed) and then using enough 
chordwise quadrature points to ensure an acceptably accurate result. An 
obvious drawback to this approach is that if a smaller value of J were used, 
more chordwise points would be used than necessary. Another, and more impor- 
tant, drawback is that actually determining the number of chordwise points to 
use requires considerable experimentation even for the regular pressure modes. 
(An example of this is provided in appendix B, figure 1 of ref. 11.) Since the 
flap- pressure mode contains additional parameters, such experimentation was not 
considered practical. Another method for determining the number of chordwise 
integration points was needed. 


T 


The method used is as follows: The chordwise integration is performed 

repeatedly with an increasing number of points and the difference obtained with 
each increase is compared until it is less than some prescribed value, S. The 
difference is approximately proportional to the actual error (when |n-n’| is 
very small this is not always true and additional considerations are necessary). 
Consequently, the actual error can be controlled by specifying 6. This proce- 
dure would be inefficient except that it is possible to nest the quadrature 
rules so that the number of arithmetic computations performed is only slightly 
greater than the number that would have been performed if the computation had 
been done only with the final number of integration points. Consider the 
chordwise integrals, which can be put in the following form; 

n, Tij) = va-x'2 h(x’, n, rij)dx’ (15) 

The above integral is evaluated by the Gaussian quadrature rule having 
/l-x' ^ as its weight function on the interval (-1, +l): 


n, ri^) - sin2 ^ h(cos n, nj 


m=i ,M 


'J' 


(16) 


Equation (l6) shows that if Gf^ resiilts from using points, then 


*^^n+l 2 *^^n M + 1 

n+1 


2 mrr , mTi 

sin"" 77 rr-h(cos 7 ; , n, n, 


M ^ + 1 

n+1 


M 


(IT) 


n+1 


if 


M = 2M +1 
n+1 n 

Thus, for example, one can obtain the result of using 15 integration points by 
computing function values at only 8 points and using also the result obtained 
with 7 points. One would then have the results of using 15 and 7 points (and 
hence information on the convergence) with only slightly more effort than would 
be involved in obtaining the answer for 15 points only. 

Now that the method whereby stability can be assured with respect to 
increasing the number of spanwise integration points (as exemplified by eq. 
(l4)) and the method for controlling the accuracy of the chordwise integration 
has been discussed, it is appropriate to consider how to minimize the total 
error for a fixed number of spanwise integration points while holding the total 
number of integration points constant. In doing this, some crude approxima- 
tions must be introduced. However, the equations representing the approxima- 
tions will be eliminated so that the results will be essentially independent of 
them. The only effect of the approximations would be to control the total num- 
ber of integration points used in obtaining the downwash values, that is, the 
procedure may not actually be fully optimized, but the results will be 
accurate. 
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(18) 


The application of equation (ll) gives 




“fj 


The right-hand side is to be minimized with respect to variations of the number 
of chordwise points used at each station, M j , subject to the constraint that 

the total niunber of integration points used be constant, that is. 


M 




M. 


(19) 




should be constant. Lagrange's multiplier method gives the following equations 

for M . ; 

J 




8M. 


+ X = 0 


( 20 ) 


Here the approximation previously discussed is introduced, namely. 


60 . 




mY 

J 


( 21 ) 


where q,r = constants > 0. 
increases. The effect of 


It is only reasonable that 6Gf 


decrease as 

bj on the error occurs because the error obtained 
with a given number of chordwise integration points increases with decreasing 
distance of the control point from the integration point due to the nature of 

K in equation (2) and figure 5. The 



Figure 5.- Illustration of the discon- 
tinuity of the function K (eqs. (2) 
and (Alb)). 


quadrature weights bj also increase 
in magnitude as the integration point 
approaches the control point. Thus, 
|bj I is a measure of this distance 
effect (even though it may be impre- 
cise) and this explains its appear- 
ance in equation (2l). 

The use of equation (2l) in (20) 
then gives 

K I <1+1 


M. 


I r+1 


~ const 


using equation (2l) again gives 




-- const 


( 22 ) 
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The values of q and r have not been separately studied, only their combina- 
tion as given in equation (22). A value of 1/2 (values of 0, 1/2, and 1 were 
used) gives the most accurate answers for a given total number of integration 
points. 

At this point it would be beneficial to summarize the concepts presented 
thus far: 

(1) The chordwise integral is done with an increasing number of chordwise 
points using nested quadrature rules for efficiency until the resulting differ- 
ences are less than a prescribed value, 6. 

(2) To use the fewest total number of integration points for a given num- 
ber spanwise, 6 is taken to be inversely proportional to the square root of 
the absolute value of the associated quadrature weight. 

(3) For convergence with respect to the number of spanwise integration 
points, 6 is taken to be inversely proportional to (j+l)^, thus 


lb^p/2(j+l)2 

Note that the error involved in applying this procedure is highly indepen- 
dent of any specific planform. The effect of planform changes is to change the 
number of chordwise integration points needed to actually achieve the 6 given 
by equation (23) and to change the number of spanwise integration points needed 
to make equation (ll) an accurate approximation of equation (lO). It is there- 
fore possible to prescribe 6 q for all plan forms , leaving J as the only 
parameter that must be varied to account for planform changes. Furthermore, 
since it is possible to nest the spanwise quadrature rules in a manner similar 
to that used in the chordwise integration, downwash modes computed with various 
values of J can be determined simultaneously. This allows a rapid determina- 
tion of whether the largest value of J used was sufficiently large. Thvis, 
multiple computer runs to assess the effect of J are lisually \annecessary. 


ILLUSTRATIVE COMPUTATIONS AND COMPARISONS 


This method for treating partial span flaps has been used with an extended 
version of Wagner’s lifting surface computer program (ref. 6) to compute pres- 
sures, span loading, and control surface derivatives on various configurations. 
The primary extensions made to this program (hereafter referred to as WAGNO) 
were: (l) to allow use of an arbitrary number of chordwise control points (Dr. 

Wagner’s extension) and (2) to allow use of an arbitrary downwash distribution. 
In this program, a total of N lateral control lines are placed on the wing 
leading and trailing edges and at the chordwise locations (given by eq. (A12)) 
and M strearawise control lines are distributed at the Multhopp locations 
(given by rijjj of eq. (5)). The intersections of these lines are the control 
points, where the boundary conditions are enforced. The number of control 
points is N X M, but symmetry or antisymmetry is used to reduce the actual 
number of unknowns to about half that value. Before discussing any final 
results, however, the character of the residual downwash should be examined. 
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Aspect ratio = 4.000 
Taper ratio = .600 



Figure 6.- Relative dimensions of the 
wing tested in reference 12. 


-41 I I I 1 I 

’0 .2 .4 .6 .8 1.0 

x/c 

Figure 7.- Chordwise distribution of 
given and residual downwash at 
n = 0,4457. 




Figure 6 shows the planform of a wing (ref, 12) for which detailed compu- 
tations of the residual downwash were performed, A chordwise distribution of 
the given downwash, a, and the residual downwash, a-a^, are shown in figure T. 
Note that a-a^ has a discontinuous derivative at the hinge line which could 
probably be eliminated, but the distribution is felt to be smooth enough for 
most purposes. Within 1-percent c of the hinge line, cannot be calculated 
accurately with the present computer program because of a limitation on the 
maximum number of spanwise integration points (5^3), The number of chordwise 
control stations can be chosen to avoid the necessity of computing in this 
region, however. As shown in reference 13, there are similar regions near, but 
excluding, the wing edges. The latter regions are small enough that they have 
caused no problem for N < 9 (the maximum number of chordwise control stations 
allowed by WAGNG) , Figure 7 also shows the chordwise control stations used for 
several values of N and corresponding values of a-a^ given to the lifting 
surface program. 


Spanwise distributions of 
tmuities in the derivative of 


a and a-a^ 
Ot— 01 ^ at the 


are shown in figure 8, The disc on- 
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eliminated also but, again, it is felt 
that this distribution is smooth enough 
for practical purposes. This figure 
also shows the singularity in downvash 
which will generally occur at any sta- 
tion for which the local sweep angle 
changes abruptly (in this case, at 
n = O). This singularity occurs for 
the regular pressure modes (eq. (h)) 
also and the only way it has been elim- 
inated is by introducing artificial 
rounding to the planform. However, if 
control points are not placed at or 
very near this point, no serious diffi- 
culties seem to arise and, apparently, 
the use of rounding is iinnecessary . 

This is the approach that has been used 
here. Figure 8 also shows the values 
of a-oif actually used in WAGNG for 
the fourth (N = 5 ) chordwise control 
station (x/c = 0.8536) at the 8 span- 
wise control stations used for M = l 6 . 

Chordwise residual pressure distri- 
butions calculated by WAGNG from a-af 
are shown in figure 9 for M = l 6 and 
N = 5, 6 , and 7. The differences are 
minor and support the previous state- 
ment that eliminating the discontinuous 
derivative of a-a^ is not necessary 
for practical purposes. 


Calculated and experimental 
results are compared in figure 10, The 
constant pressure panel calculation was 
performed by Ralph Carmichael of Ames 
Research Center using a version of the 
Ames Wing-Body Computer Program (ref. 
ik) . Two-hundred panels of unknown 
strength were used on the semispan. WAGNG with N = T and M = l6 used 56 
unknowns (although a comparable result was obtained with N = 5 and thus 40 
unknowns). The theoretical results agree reasonably well with experiment 
except aft of the hinge line, where the influence of viscosity is strong. 


Figure 8,- Spanwise distribution of 
given and residual downwashes at 
the chordwise station x/c = 0.8536. 


Figure 11 shows a predicted section lift coefficient distribution for a 
10° flap deflection on an AR = 8.56 swept wing. It is compared with experi- 
mental values and values predicted by a vortex lattice method for a nearly 
identical wing (ref, 15). The difference between the theoretical methods is 
small, but both differ considerably from the experimental values. Reference 15 
suggests that a possible reason for this may be separated flow since the data 
were taken at a wing angle of attack of 10 ° and thus the flap was actually at 
20 ° incidence. 
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ACn-AC 


Cl 



Figure 9.- Comparison o£ residual 

chordwise pressure distributions Figure 10.- Comparison o£ experimental 

calculated by the Wagner computer and calculated li£ting pressure 

program £rom the residual downwash distributions on a swept wing with 
a - a^; n = 0.4457. partial span £lap; M = 0.6, n = 0.46. 



0 .2 .4 .6 ,8 1.0 


y 

b/2 

Figure 11.- Comparison o£ calculated and experimental section li£t coe££icient 
distributions £or a wing with a partial span £lap. 
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Figure 12.- Comparison of predicted 
and experimental AC for a 60® 

delta wing. ^ Figure 13.- Wing tested in reference 17. 

A third comparison is given in figure 12, where ACl induced by flap 
deflection on a 60^ delta wing is shown as a fimction of angle of attack. A 
sufficiently accurate answer was obtained with N = 7, with the predicted slope 
only 3 percent above the experimental value (ref. l6). This close agreement, 
however, may not be the result of a good theory, but more likely because the 
theory ignores the presence of the fuselage', which acts somewhat as a reflec- 
tion plane and thus would tend to increase • This phenomenon is also 

apparent in the next comparison. 

Figure 13 shows the planform of this next example (ref. 17)* The Mach 
number variations of the predicted and experimental control surface derivatives, 

Ik 



Cl^, Chj, and , are shown in figures through l6, respectively (the theory 

is invalid at M = 0.9, however, since this is above the critical Mach number). 
The term is actually underpredicted by 5 percent (8 x 24 control points 

at M = 0.6), which, as mentioned above, results from not accoionting for the 
interference of the fuselage. The term Cjj^ is overpredicted by 15 percent 
(M = 0.6). Had the body interference been accounted for, the discrepancy might 
have been greater; this discrepancy is large (and expected) since depends 

only on the lifting pressure in the region most adversely affected by viscosity. 
The rolling moment derivative, is underpredicted by 6.5 percent (M = 0.6). 


The comparisons presented here indicate that the present method can pre- 
dict with reasonable accuracy features of the flow about wings with deflected 
control surfaces. However, this linearized theory will fail when there is 
separated flow, and this happens more readily when a control surface is 
deflected. Consequently, the method should be used with discrimination. The 
comparisons also indicate that finite element methods can predict flow features 
for this type of configuration. However, a method of the present type (i.e., 
employing continuous loading functions) can be made more economical since far 
fewer unknowns are needed. 
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Figure 14.- Variation of control 
surface lift effectiveness with 
Mach number (symmetric deflection 
of two surfaces) . 


Figure 15.- Variation of the 
hinge moment coefficient with 
Mach number (one surface only) . 
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Mach number 


Figure 16.- Variation of control surface rolling moment effectiveness with 
Mach number (antisymmetric deflection of two surfaces) . 

CONCLUDING REMARKS 


The governing integral equation of linearized, steady lifting surface 
theory was examined and a standard method for solving it, based on infinite 
aspect ratio theory, was briefly reviewed* It was pointed out that to treat 
the partial span flap, one should not proceed similarly from the corresponding 
infinite aspect ratio theory but rather use a three-dimensional local solution 
(ref. l) because only the latter can reliably represent the pressure distribu- 
tion near the control surface corner. This local solution was incorporated 
with an existing lifting surface program to obtain a complete solution to the 
problem. 

Next a basic procedure for performing the numerical quadratures required 
was derived. The results were in no way related to the fact that a control 
surface pressure mode was being used, so they actually bear a direct signifi- 
cance to lifting surface theory in general. It was determined how instabili- 
ties in the spanwise integration occur if the chordwise and spanwise integra- 
tion procedures are not coupled. The manner in which the integration procedure 
should be related required the chordwise integration to be performed with a 
specified accuracy, depending on the spanwise integration and control stations. 
It was indicated that the usual procedures are necessarily inefficient in 
attaining the required accuracy, A new method was introduced which is more 


l6 



efficient because it does not compute to a greater accuracy than needed and, 
furthermore, because it uses the concept of ^nested” q_uadrature rules. 

Following the discussion of the q^uadrature procedures, illustrative compu- 
tations and comparisons with experiments were discussed. The nature of the 
control surface downwash mode was shown to be such that the discontinuity in 
the given downwash could be properly eliminated, but discontinuous derivatives 
remained. Calculated examples showed, however, that these irregularities are 
sufficiently weak so as not to constitute a real problem. The method presented 
was demonstrated to be an effective, accurate solution to the integral equation. 
However, in some cases there were rather large discrepancies with experiments. 
This was not unexpected. The principal causes were viscosity and the neglect 
of body and thickness interference effects. Also, another possible source of 
discrepancy is the gap between the control surface side edges and the wing due 
to the control surface’s deflection. This gap is not accounted for in the 
theory. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 9^+035, July 26, 1972 
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APPENDIX A 


FORM OF THE FLAP PRESSURE MODE AND CHORDWISE INTEGRATION 


The flap downwash mode is given by 




te ACt 


where 


he 


K[C-C* ,3(n-n' )]dC dn' 


(Ala) 


K[c-5», 3(n-n’)] = 1 + 


irlL 


/(5-r)2 + 62(r,_n')2 


(Alb) 


This can he put in a more convenient form hy the use of the following 
transformations : 


’te' 


then 


where 


’te' 


* + l 



(A2a) 


(A2b) 

^ Tf K(x-x* , y)dx* dn’ 

(A3) 


h3 ( t \ 

Expressions for ACp were given hy Landahl and extended hy Ashley (ref. 
l8) to the case of a swept flap. The expression given helow (for unit flap 
deflection) has heen adapted from these references with a correction made to 
properly account for the flap edges: 


where 


-2 cos A 


AC = 


^ tt/1-M2 cos^X^ 


F(x’ ,n’ )i^n 


/(x’-x^)2 + - yjL 


_/(x*-X^)2 + y2 - y^ 


^ /I - M2 cos2X 

= ■Jirrn’ — ^ 


(AUa) 


(AUb) 


•u /I - m2 cos2x 
b c 


c(n ’ ) cos 


(r)*-ri2) 


(aUc) 
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The function F(x’,n*) must attain a value of 1.0 on the flap hinge line and 
must introduce the correct edge behavior on the periphery of the planform (ref. 
l), but is otherwise arbitrary. It has been chosen to be the following: 


F(x»,n’) = Co(n’) 


a - X 


i2 


1 - x„ 


where 



2-1 


1^2 


1 1 in' |-m2i 


l-|ri2 


(AUd) 


n’ < ri2 


n* I > |ri2 


(Ake) 


Combining equations (A3) and (aU) yields 


where 


and 


cos X 


a^(5,n) = 

4ir^/l-M2 cos2X 


4 

c *Li 


« C,(n') ^ 


(n'-n)" 


n, n*)dn» (A5a) 


+1 


H^(C, n, n*) = / ^ /l-x‘2 g(x*)<3x' 


g(x» ) = — i — )i,n 
c 


v'(x'-x^)2 + y^ _ y^ 


/(x’-X )2 + y2 _ 


'2J 


K(x-x’ , y) 


(A5b) 


(A5c) 


As indicated previously, the function may have to be calculated very 

accurately and it is to be calculated by the following Gaussian quadrature 
rule: 


/l-x*2 g(x')dx' = 


where 


j+i 


(A6a) 


(A6b) 


For equation (A6a) to be accurate, the function g(x*) should be able to 
be expressed accurately as a polynomial of degree less than 2J. Accordingly, 
the logarithmic singularity of g(x* ) at x^ when Yi and Y 2 are of opposite 
sign or when Yi or Y 2 is zero mxist be eliminated or at least softened. Simi- 
larly, the discontinuity that K has at x' = x when j’’ = 0 must be softened. 
The former has been accomplished with the following fxmction: 
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which is not singular anywhere, hut has singular derivatives at x’ = ±1 and 
x’ = , hut these are much less serious. The influence function is then given 

by 




When y = 0 and |x| ^1, there is a discontinuity within the integration 
interval (fig. 5). The previous procedure leads to the following expression 
for : 

K(x-x ,0) ^ 

/l^x’ 2 g(x’ )dx + {same as in (a8) . . . } (AS 

-1 

The numerical integration is accomplished hy substituting 


g(x* ) = 


(x' ) - g(x) - (x'-x)g'(x) x' < 


x' > X 


This function is suitably smooth to be integrated by the Gaussian 
quadrature method, resulting in: 


20 


(All) 


H^,(C,n,n*) =1^ ^ (l-x )g(x ) + g(x)[x/l-x'2 + sin“^x + 

f J+1 J J 2 


-g* (x) 


2+x^ 


/l-x2 + x(sin”^x + 


+ last term in (A9) 


In the above form^a, note that the summation need not be performed for Xj > x 
and that g(x) = g'(x) = 0 if x < x^. 

A question deserving consideration is whether there is a best choice for 
the number of quadrat\ire points, J. Consider the function K(x-x* ,y) (see fig. 
5). For small values of y, K(x-x* ,y) is nearly discontinuous at x-x* =0, 
where it attains an average value. It seems obvious that an integration point 
at X* = X would be desirable since sane of the errors could then be made to 
cancel. The lifting surface program with which this method is used has (for 
small Y) the control points at x = ±1 and 

X = cos p = 1,2,. . ., N-2 (A12) 


From equation (A6b), then, appropriate values for J would be 

J = (N+l)^ - 1 


(A13) 


where k is any integer • 

A final point to mention is that the procedures given for determining 
Hf(5jn,n’) may also be used to determine the flap pressure mode contribution to 
the section coefficients of lift and pitching moment about the quarter-chord 
merely by replacing K with 1 and with x*-0,5 5 respectively. For this compu- 
tation, however, there is no optimum value for J. 


21 


APPENDIX B 


CORRECTIONS FOR THE LOGARITHMIC SINGULARITY AND OTHER SINGULARITIES 

IN THE SPANWISE INTEGRAL 


The spanwise integral for the flap downwash mode contains a factor of the 


form 


H(C 


,n,n') = 


r+1 

1 h(x' ,n' ) 

+ 

J 

J-1 

L /(x-x’ )2 + y2J 


dx' 


(Bl) 


where h(x’,n’) is a function that behaves as /l-x' 2 near the end points and 
that may have a logarithmic singularity. Assume, however, that in the vicinity 
of X* = X it is a regular function. By developing h(x*,n') in a Taylor 
series it can be shown that H(^,ri,ri') has a logarithmically singular second 
derivative. Only the contribution in the vicinity of x* = x need be 


considered: 


x+6 

( 


H(c,n,n’ ) = I h(x' ,n* ) 


x-x' 


x-iS 

x+6 

X 


/ ( x-x * ) ^ + y ^ 


9h 


dx' + a regular part 


(h(x,n') + (x'-x) (x,n') + . . . 


x-x' 


^-6 

a 

■+6 

-6 


/(x-x' )2 + y2 


dx' 


+ a regular part 
9h 


i 


9x 


7 (x,n’ )x + . . . 


-X 


/x2+y2 


dx + a regular part 


o f I \ 

= -2 (x,n’) 


^ /x2+y2 _ ^ Jln(x+ v^+y^) 


+ a regular part + lower-order singularities 

= - (x,n’ )y^ £n|Y| + a regular part 

+ lower-order singularities 

Vk ft 

' ^ I ’1“'^ * I S' i*eg^ar part 


9x 


r (x,n' ) 


+ lower-order singularities 


(B2) 


The factor (n-n')^ is cancelled by its reciprocal in the spanwise inte- 
gral, leaving a logarithmic singularity. This singularity also exists for the 
regular downwash modes and was recognized by Multhopp (ref. 2), who devised a 
method to correct for it. Subsequently, Mangier and Spencer (ref. 19) intro- 
duced a better method and Garner (ref. 20) illustrated the importance of this 
correction. This method of correcting for the logarithmic singularity has been 
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used for the flap pressure mode and is reproduced here. Only the most impor- 
tant terms are retained (i.e., ( 3h/9x' ) (x,n’ ) and c(ri') are replaced by 
(9h/9x)(x,n) and c(n)). The proper edge behavior must be int roduced into the 
singular part; accordingly, it is multiplied by the factor /l-n *2//l-n2 . The 
singularity is then added and subtracted from the integrand, giving the follow- 
ing expression for the flap downwash mode: 


a^(5,n) 


cos X 


1+tt2/1-M2 cos2x. 


f 


C 2 (n') 


c(n’ ) 


A-r )2 


+C 2 (n )32 (x,n) 

C2(n)32b 


(n’-n )2 

Jlnlri-n* 




dn' 


A-n2 c(n) 


1^ (x,n) / ^ A-n'2 j),n|ri-n' Idn’y 


(B3) 


The bracketed term will not contain a logarithmic singularity and is inte- 
grated numerically by M\ilthopp's integration formula. The integral in the last 
term has been done analytically in reference 20 : 


/ ^ /l-n ’ 2 5,nln'-n|dn = ■^[ 2 n^ - l - tn 4] 
The expression for h(x',n') is 

/(x*-x^)2 + y^ - y^~^ 
/(x'-x^)2 + y^ - y 2 


h(: 




l-x'^ 

l-x2 

c 


tn 


(b4) 


(B5) 


which, when combined with equations (B3) and (bU), will yield the final form of 
the equation for the flap mode of pressure for an interior control point. For 
a control point on a wing leading edge (x’ = -l) or trailing edge (x’ = +l) , 
however, 9h/3x’ becomes infinite, violating the assumption that h(x’,nO can 
he developed in a Taylor series. 


The appropriate approach in either instance is to retain the square-root 
behavior in the chordwise integral and to develop the remainder of the function 
h(x’,n’) in a Taylor series about the leading- or trailing-edge point. This 
approach was used by Jordan (ref. 2l) for regular pressure modes. When applied 
to the flap pressure modes , the following expression has been derived for the 
lowest order irregular term in the flap influence function for the leading 
edge: 




16 . 

*^(l+Xc)2+yf - y^ 


bg 

3A-x2 
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/(l+X^)2+y2 _ 


2 c(n)cos 


3/2 


(B 6 a) 
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where 




V = I <1 ■* +15) 


'£e 


= tan' 


.-1 


-r tan 
e ie 


(b6i3) 

(b6c) 

(B6d) 


where is the leading-edge sweep angle at the control point defined so 

that -it/2 < < +it/2 and X^^ > 0 on that side of the control point for which 

the edge sweeps forward. The leading-edge sweep angle is of the opposite sign 
on the other side (assuming no kink in the leading edge). The corresponding 
formula for the trailing edge is 




where 
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3/l-x2 
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’^(i-Xc)^+yf - yj 


/(l-X )2+jr| - 72 


he 

2c (n) COS (f) 


te 


3/2 


”^te = I (1 - 


*^te 


= tan“^ 


— tan X. 

e te 


Fi(m^^)|n'-n|3/2 
(BTa) 
(B7b) 
(BTc) 


The term X^g is defined as for the leading edge. Note that in either 
case the parameter, m, of the elliptic integrals changes on passing from one 
side of the control point to the other (more precisely, if m"*" is the value for 
n’ > n, then 1-m"*" is the value for n’ < n if there is no kink at n’ = n). 

As in the case o f int erior points, this irregularity is, after introducing the 
factor A-n to account for the edge behavior, simultaneously added 

and subtracted from the integrand, leaving a part to integrate analytically and 
a smoother part to integrate numerically. At the leading edge, equation (B6a) 
is actually the leading term, so the spanwise integral does not have the 
second-order pole, but rather one of order -1/2, After this is eliminated by 
adding and subtracting the irregular part (eq, (B6a)), the leading term of the 
integrand in the vicinity of n* = n would be a constant followed by terms of 
order 1/2, 1, etc. Multhopp^s integration formula may still be conveniently 
used for this integration. The integral for the leading edge is thus 
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(equation continued on the following page) 
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(b 8 ) 


The bracketed term is not singular and is computed using Multhopp’s formula. 
The last term, computed analytically, is also expressed in terms of elliptic 
integrals : 
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Therefore, for the leading edge, 
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For the trailing edge, 



(Bll) 

It is important to note that values for , yi , 72 , 4 >„ > and (|i appear- 

jce "C e 

ing in the above equations are to be evaluated at the control point, n. An 
appropriate question might be: ”How important are these corrections?" The 
answer is supplied by figures lT(a) through (e), which show downwash values cal- 
culated for various points on a sample wing for various numbers of spanwise 
integration points, J, These graphs clearly show that it is important to use 
these corrections, at least when Multhopp’s integration formula is used. For 
other integration procedures (i.e., the zonal scheme of Watkins et al. (ref. 
22 )), this method of correction could also be applied. Although it is not 
known how important the corrections would be, it seems, in view of how well 
they worked for Multhopp’s method, that they could be significajitly worthwhile. 
Figure 17 also shows how close a^(5,n) can be predicted. An appropriate num- 
ber to compare these values with is 1, since the flap pressure mode produces a 
discontinuity of magnitude 1. For J = 191, the error is less than 0.3 percent 
for all points, except for the one which is only l.U percent c behind the 
hinge line. The error for this point is approximately 6 percent. When comput- 
ing ot^(^,n) for the lifting surface program, however, N can be chosen so that 
all or most of the collocation points are not so near the hinge line. The 
spanwise control point positions are chosen so that they are not near the side 
edges of the control surface. 
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Figure 17.- Continued. 
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